Synchronization in clustered random networks 
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In this paper we study synchronization of random clustered networks consisting of Kuramoto 
oscillators. More specifically, by developing a mean-field analysis, we find that the presence of 
cycles of order three does not play an important role on network synchronization, showing that the 
synchronization of random clustered networks can be described by tree-based theories, even for high 
values of clustering. In order to support our findings, we provide numerical simulations considering 
clustered and non-clustered networks, which are in good agreement with our theoretical results. 
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I. INTRODUCTION 



Synchronization processes have attracted the interest 
of scientists for centuries and is in the focus of intense 
research today This collective phenomena has been 
observed in biological, chemical, physical, and social sys- 
tems 0, [1| • Many works have verified that the dynamics 
of synchronization depends on the connectivity pattern 
of networks For instance, when the natural frequency 
distributions are unimodal and even, the critical coupling 
depends on the ratio between the first and second statis- 
tical moments of the degree distribution ■ In addi- 
tion, for networks in which there is a positive correlation 
between the network structure and dynamics, the criti- 
cal coupling has an inverse dependence on the network 
average degre 

However most of the analytical results of networks con- 
sisting of Kuramoto oscillators have been obtained for un- 
correlated networks generated through the configuration 
model [J-Q , which generates networks with arbitrary de- 
gree distributions by randomly connecting the nodes ac- 
cording to a specified degree sequence. One of the main 
properties of this model is that in the thermodynamic 
limit; i.e., N — > oo, the probability of occurrence of cy- 
cles of order three tends to zero. Such probability can be 
quantified through the clustering coefficient C, defined 



as 



C 



3 X (number of triangules in the network) SiV, 



A 



number of connected triples 
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In addition to the configuration model, the cluster- 
ing coefficient also vanishes for Erdos-Renyi (ER) and 
Barabasi- Albert (BA) networks when iV — > oo. There- 
fore, most of the theoretical results concerning the Ku- 
ramoto model are derived for networks that have locally 
tree-like structures, i.e., C — > 0. However, most real- 
world networks have clustered topologies Thus, an 
analysis of the Kuramoto model on networks with clus- 
tering is necessary to model real-world synchronization 
with more accuracy. 

The objective of the current work is to study the dy- 
namics of synchronization on clustered random networks 
and perform a comparison with non-clustered ones. More 
specifically, we develop a mean-field theory for the con- 
figuration model proposed independently by Newman Q 
and Miller [10(], which generates networks with C > 
even in the limit of large networks. Such analysis is com- 
pared with the mean-field theory developed for locally 
tree-like networks. Our results show that the mean- 
field theory for networks with low values of triangles 
can be applied with certain accuracy on clustered net- 
works, indicating that the presence of cycles of order 
three does not influence the network synchronization. 
This result is in agreement with previous works 
which observed that the clustering coefficient does not 
play an important role in other dynamical process, such 
as bond-percolation, rumor and epidemic spreading, pro- 
vided that the networks have low values of the average 
shortest path length. 
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In Sec. II we briefly describe the configuration model 
proposed in and [l^- In Sec. Ill we derive a suffi- 
cient condition for synchronization through mean-field 
approximation and in Sec. IV we compare numerical and 
theoretical results and give our conclusions. 



where Wi is the natural frequency of the node i , A is the 
coupling strength and Aij are the elements of the adja- 
cency matrix A, where Aij = 1 if the nodes i and j are 
connected while Aij ~ 0, otherwise. The synchronization 
can be quantified through the order parameter 



II. RANDOM CLUSTERED NETWORKS 

In the standard configuration model, the network is 
generated through the degree sequence {ki}, connecting 
the "stubs" at random. The process to generate random 
clustered networks is quite similar. Let Si and ti be the 
number of single edges and the number of triangles at- 
tached to the node i, respectively. Given a network the 
sequence {si,ti} is possible to connect the "stubs" in or- 
der to generate single edges and also to connect nodes in 
order to obtain complete triangles. Hence, it is conve- 
nient to define the joint degree distribution Ps.t^ of the 
network, which is the fraction of vertices connected to 
s single edges and triangles. Therefore, the conven- 
tional degree of each node is given by fc = s -I- 21^, since 
each triangle contributes with two to the degree. Also, 
it is possible to relate the joint degree distribution Pst& 
with the conventional degree distribution pk through 



Pk 



E 

s.t=0 



Pstzs0k,s+2t&, 



(2) 



where Sij is the Kronecker delta. 

With the joint degree distribution Ps,tA ^^'^ the degree 
distribution pk , we can calculate the clustering coefficient 
for random networks. The number of triangles in the 
network is given by STVa = X!st ^Pst ^^'^ the number of 



connected triples 7V3 = Nj^k 
the clustering coefficient is 



)pk- Thus, using Eq. [U 



C 



Est ^P«tA 

Ek 



(3) 



Note that the of factors N cancel, letting C > in the 
limit of large networks, i.e., N — > 00. 



III. SYNCHRONIZATION ON CLUSTERED 
NETWORKS 

The Kuramoto model consists of a set of N oscillators 
coupled by the sine of their phase differences 0, [l^ . The 
state of each oscillator is characterized by its phase 9i {t) 
i = 1, . . . , N . Considering a complex network where each 
node is a Kuramoto oscillator, the equations of motion 
are given by 



dt 



N 



sm 



-e,), i^i,...,N, (4) 



re 



1 



N 



N ^ 



(5) 



where ipit) is the average phase of the system. The 
coherence parameter is bounded as < r < 1, where 
r = 1 represents the fully synchronized state and r = 
is the incoherent solution. In the fully connected graph 
(Aij = 1 Vi, J and i ^ j), the order parameter r as a func- 
tion of A displays a second-order phase transition char- 
acterized by the critical coupling Ac = 2/(7rg(w)) (sl. [l5|. 
where g{Lo) is the distribution of the natural frequencies 
and u) is the average frequency. In random networks, the 
critical coupling Ac of such phase transition is rescaled 
by the ratio (fc) / (fc^) i.e.. 



Ac = 



(fc) 



iTg{uj) (fc2) ' 



(6) 



where (fc") is the n-th moment of the degree distribution 
Pk of the network. 



A. Mean-field theory 

The critical coupling strength necessary for the onset 
of synchronization in Eq. [5] was firstly obtained by Ichi- 
nomiya through a mean-field analysis for the standard 
configuration model. The mean-field analysis has the ad- 
vantage of allowing an analytical treatment. In order 
to extend the mean-field treatment to the configuration 
model for random clustered networks, we approximate 
the Eqs.|4]into the following equation 



dt 



CO, + Xy hP{k'\h) sin {9k' - e,) , (7) 



fe 



where ki is the degree of the node i and P(fc'|fc) is the 
probability that an edge emitted by a node with degree fc 
is connected to a node with fc'. For uncorrelated networks 
P(fc'|fc) is given by 



P(fc'|fc") = 



k'pk' 
{k) 



Substituting Eq. |S] in Eq. H] we have: 

de,{t) , Afc, 



dt 



(fc) 



Pk' sm ( 



(8) 



(9) 



For the configuration model of clustered random net- 
works, Pk is defined by Eq. [51 Therefore, substituting in 
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FIG. 1: Synchronization diagrams for networks with double 
Poisson joint distribution Pst^ (Eq. [23} • The dots are ob- 
tained calculating the equations of motion (Eq. [4} until the 
system reaches the stationary state for each value of coupling 
A. The order parameter r is then calculated with Eq. [5] Each 
point is an average over 10 network realizations. Solid lines 
correspond to the theoretical prediction from Eq. 1211 



Eq. [H] and noting that ki = Si + 2ti we get 



dt 



UJ,- 



^ (s'+24)p,,t^ sin( 



(10) 

where (fc) = (s) + 2 (<a)- 

For an analytic treatment it is convenient to use the 
continuum limit of Eq. [TOl For this purpose, let us define 
the density of the nodes with phase 9 at time t, for a 
given w, with s singles edges and t/^ triangles, denoted 
by p(s,t/:^,u;;9,t). This density is normalized as 



27r 



(11) 



Therefore, Eq. (TU] in the continuum limit is given by 



de{t) 



A(s + 2tA) 



dt (k) 
xps't' sin {6 



ds' I dt'^ I de'{s' + 2t'^, 
(12) 



The order parameter can also be redefined in order to 
account the connectivity pattern of a random network as 

re'^^''^ = ^ JdLjJdsJ dt^ J d9is + 2tA)Pst^ 

xp{s,tA,uj;0,t)e''^- (13) 

Considering Eq.[T31 it allows to rewrite Eq.[Tn]in terms 
of the order parameter, resulting in 

do 

— =u; + \{s + 2t/^,)sm{'^lJ-e). (14) 
dt 

The density p{s,t/x,Lu; 6,t) will obey the following conti- 
nuity equation 

dp d 



dt 86 



{vep{s,ti^,uj]9,t)} = 0, 



(15) 
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FIG. 2: Synchronization diagram calculated as in Fig. [T] for 
networks with double power-law joint distribution Pst^ , with 
7s = 7t = 3 in Eq. 1241 Each point is an average over 10 
network realizations. Solid lines correspond to the theoretical 
prediction from Eq. 1211 



which for the stationary states {dp/dt = 0) has the solu- 
tions 



p{s,tA,Lo; 



— arcsm 

\uj — Xkr sin 6 \ 



" (s + 2tA) 

otherwise. 



< Ar, 



(16) 

These solutions correspond to those oscillators that are 
entrained by the mean field and those non-entrained, re- 
spectively. Thus, separating each contribution to the or- 
der parameter we yield 



{k)r 



ds I dt/xd9 



r<\r 



X Pst^is + 2t/^^)g{uj)p{s,tA;0)e''' . 



>Xr 



(17) 



The part of the non-entrained oscillator is given by 



d9g{uj)e' 



duj 



-A(s+2tA)r 



+2tA)r (w- A(s + 2tA)rsin^ 



duj- 



(— a; + A(s -I- 2tA)r sin t 



= 0,(18) 



since the integral over is 0. Thus the only contribution 
for coherence the parameter r is due to the synchronous 
oscillators, which is accounted in Eq. [iTl 



(fc) r = ds dt/s I duj \ dO{s + 2iA)g(w) 



xpst^ exp 



I arcsm 



{s + 2tA)Xr 



(19) 
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From the real part we get: 

{k)r = ds dt/s I duj I de{s + 2tA)g{uj) 



(s + 2tA)Ar 



(20) 



Considering the following change of variable w' = Ld/{s + 
2t/\)\r and considering g{u}) = (•\/27r)~^e~" we ob- 
tain the following implicit equation for the coherence pa- 
rameter r 



dsdt 



A 



(21) 
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FIG. 3: Synchronization diagrams calculated as in the 
other figures for random networks with degree distribution 
Pk = e^^'^^ /k\ and joint degree distribution generated 
by Eq. 1251 Each point is an average over 10 network real- 
izations. Solid lines correspond to the theoretical prediction 
from Eq. [m 
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FIG. 4: Synchronization diagrams calculated as in the other 
figures for random networks with degree distribution pk oc 
k^~' , with 7 = 3 and joint degree distribution generated by 
Eq. 1251 Each point is an average over 10 network realiza- 
tions. Solid lines correspond to the theoretical prediction from 
Eq.im 



where Iq and Ii are the modified Bessel functions of first 
kind. Thus, tending — > 0, we obtain the critical cou- 
pling Ac for the onset of synchronization 



Ac 



(fc) 



(s + 2tA)^pst^dsdtA 



(22) 



The mean field result obtained for the critical coupling 
Ac for clustered netwokrs, Eq. [22l is similar to the one 
for non-clustered networks (Eq. i.e., in both cases 
the critical coupling is proportional to the ratio of the 
moments of the degree distribution. Note that in the 
absence of triangles (^a = 0) we recover the result Ac = 
\/S/tt (k) I (fc^), where the degrees are just due to single 
edges, k = s. 



IV. NUMERICAL SIMULATIONS 

In this section we give some numerical simulations of 
clustered and non-clustered networks and will compare 



them with the theoretical result of Eq. [2TJ All simu- 
lations consider networks which are constructed through 
the configuration model presented in Sec. II and the distri- 
butions of single edges and triangles independently. How- 
ever, it is also possible consider correlated distributions 
as well. 

Let us study first networks with the following joint 
distribution of single edges and triangles 



Pst^ 



-{tA> 



A 



(23) 



In order to analyse systematically the dependence of the 
order parameter on the presence of triangles in the net- 
work, we kept the average degree (fc) = (s) +2 (^a) fixed 
and varied the (s) and {t/\), calculating the order param- 
eter r as a function of the critical coupling A. Fig. [T] shows 
the synchronization diagram for networks with double 
Poisson degree distributions (Eq. [23|) with average de- 
gree (fc) 20. It is interesting to note that networks with 
higher values of (t/s) have the same critical coupling for 
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the onset of synchronization. 

We have also considered networks with joint distribu- 
tion consisting of a double power-law distribution 



'■A ' 



(24) 



where = It = 1 for the sake of simplicity. Fig.[5]shows 
order parameter r as a function of A considering 7 = 3. 
As we can see, the same behavior is observed as in Fig.[TJ 
the presence of clustering in the network does not affect 
the network synchronization. The non-zero values of the 
order parameter r for small values of the coupling A in 
Fig. [T] and [2] are due to finite-size effects 

It is also possible to construct the joint distribution 
Pst^ from a g iven degree distribution pk through the re- 
lation [H, [Til 

Pst^ = Pk5k,s+2t^ [(1 - f)5tfl + f5tX{s+2t^)/2\\ (25) 

where < / < 1 and [-J is the floor function. Through 
Eq. [25] we can construct Pst^ keeping the degree distri- 
bution fixed with / being the fraction of nodes in the 
network attached to the maximum possible number of 
triangles t ^ \{.s + 2t/^)/2\ and (1 — /) the fraction of 
nodes which are attached to single edg es only. Substitut- 
ing Eq. [25] into Eq. [3] we obtained 



C = / 



Sfcbzfc +P2fc+l) 



(26) 



Eq. [55] establishes a linear relationship between C and /. 



i.e.. with / = we construct a network with the mini- 
mum value for the cluster coefficient and / = I a network 
with the maximum value of C for a given p^, allowing 
to study the extreme cases of the topology. Fig. [3] and 
Fig. [5] show the synchronization diagrams for networks 
with pk = e"^'""^ {k)^ /k\ and pk oc ^ respectively. 
Again, we observe a good agreement with the theoret- 
ical curve. Therefore, the clustering coefficient has no 
effect on the coherence parameter evolution r(A), com- 
paring the curves with / = and / = f . Also, finite-size 
effects are observed for small values of A. 

In summary, we have shown that the presence of cycles 
of order three does not play an important role in network 
synchronization of Kuramoto oscillators. In fact, the the- 
oretical results for non-clustered networks are highly ac- 
curate on describing the behavior of the order parameter 
r for clustered networks, even when the cluster coefficient 
C has the maximum accessible value for a given network. 
The results presented here are in agreement with previ- 
ous findings [ff|, where it was found that the presence 
of triangles in the network topology does not influence 
the performance of other dynamical processes, such as 
bond percolation, fc-core size percolations and epidemic 
spreading. 
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